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. We propose a new quantum dynamics method called the effective potential analytic continuation 

\ (EPAC) to calculate the real time quantum correlation functions at finite temperature. The method 

is based on the effective action formalism which includes the standard effective potential. The basic 
notions of the EPAC are presented for a one-dimensional double well system in comparison with 
the centroid molecular dynamics (CMD) and the exact real time quantum correlation function, ft 
is shown that both the EPAC and the CMD well reproduce the exact short time behavior, while at 
longer time their results deviate from the exact one. The CMD correlation function damps rapidly 
with time because of ensemble dephasing. The EPAC correlation function, however, can reproduce 
the long time oscillation inherent in the quantum double well systems. It is also shown that the 
EPAC correlation function can be improved toward the exact correlation function by means of the 
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O ' I. INTRODUCTION 



The imaginary time path integral Q provides the quantum statistical mechanical formalism of useful numerical 
methods for quantum many-body systems. Static observables in such systems have successfully been calculated by 



^J^- numerical methods have difficulty in accessing the dynamical properties such as the real time quantum correlation 



the path integral Monte Carlo (PIMC) or path integral molecular dynamics (PIMD) calculation y,yj. However, these 



functions, because the analytic continuation from the imaginary time to the real time U,|5|,|6| using a finite number of 
noisy imaginary time PIMC/PIMD data is so nontrivial as to be classified as an ill-posed problem. To overcome such 
difficulty, the numerical analytic continuation scheme based on the maximum entropy method has been proposed to 



be applied to various quantum dynamical problems [3, la l2| ■ As an alternative approach, one can directly evaluate the 
real time path integrals, although this suffers from the sign problem; the err or g rows exponentially with time because 
of the rapid oscillation which originates from the factor exp(iS[q(t)]/h) |2|, |lfj, |llj] . In addition to these approaches, 
another numerical method has been proposed to evaluate the eigenstates of quantum systems using the path integral 
and to construct the real time quantum correlation functions |l21 . 

The effective potential is a device widely used for the approximate calculations of static properties of quantum 
many-body systems. There are several definitions of effective potentials: effective classical potential [l|, 1 13 . Il4j . 
standard effective potential 0| , and so on [lfl \v\ fisl Q] . All of them are defined in the framework of path integral 
while the relations among them have been discussed so far (lii Efl I21I E3 . Ej] . The usefulness of the effective potentials 
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has been indicated by many applications to simple quantum mechanical systems 0, 0] , condensed phase systems 
[2^ , the quantum transition-state theory [25! Eg , and quantum field theories [2^ . It is true that static properties 
of the systems are well described in terms of the effective potentials in a simple classical analogue. For instance, 
the quantum-mechanical partition function is expressed as a classical-like partition function including the effective 
classical potential. Furthermore, both of the thermodynamic phases of quantum statistical systems and the vacuum 
structure of quantum field theories are determined by the minima of the standard effective potentials. However, it has 
been believed that such classical use of the effective potentials cannot be directly applied to the dynamical problems 
because the definitions of the effective potentials, in most cases, do not justify the use of the classical equations of 
motion on the effective potential surface. Therefore, for the calculation of time-dependent properties such as the real 
time quantum correlation functions, more careful treatment of the effective potentials is required. 

On the other hand, Cao and Voth have recently proposed the centroid molecular dynamics (CMD) method, which 
is an approximation to obtain real time quantum correlation functions from molecular dynamics on the effective 
classical potential surface [28^ . The validity of this method is ensured by the fact that the real time correlation 
function of the centroid variables is a good approximation to the exact canonical correlation function in the linear 
response theory 00. The CMD is a promising method suited for the numerical computation of the dynamics of 
many-body molecular systems such as condensed matter and molecular clusters [ilj. In fact, the validity of the CMD 
approximation has been tested for low dimensional nonlinear systems; the real time correlation functions obtained 
from the CMD are more accurate at shorter time and at lower temperature, whereas they are evidently worse at 
longer time It has also been found that the time correlation function evaluated by the CMD damps rapidly 

with time because of ensemble dephasing, which is well-known behavior in one-dimensional nonlinear classical systems 
30|. 

In the present paper, we propose a new quantum dynamics method, called the effective potential analytic continu- 
ation (EPAC) method, to calculate the real time quantum correlation functions. The EPAC method is based on the 
effective action formalism 00Q00, where imaginary time quantum correlation functions are expressed in 
terms of the effective action. This method is an approximation method which includes the standard effective potential 
defined as the leading order of the derivative expansion of the effective action. Once the standard effective potential is 
known, one can obtain the analytic form of the imaginary time quantum correlation function, and then the analytic 
continuation from the imaginary time to the real time can be readily performed. The EPAC is expected to be a 
powerful method to calculate the real time quantum correlation functions with less computational effort than the 
other numerical methods, because the standard effective potential can be easily calculated from the effective classical 
potential |22( . 

In the present work, we apply the EPAC method to a one-dimensional quantum double well system at finite 
temperature. At first, the standard effective potential is obtained from the effective classical potential calculated by 
means of the PIMD technique. And then the real time quantum correlation function is constructed by means of the 
analytic continuation procedure. The results are compared with the exact quantum statistical mechanical results to 
investigate the accuracy of the EPAC method. We also compare the EPAC results with the CMD results, to clarify the 
difference between the approximate quantum dynamics based on the standard effective potential and on the effective 
classical potential. Finally a possible improvement of the EPAC is discussed on the basis of the derivative expansion. 

In Sec. II, we summarize the definitions and fundamental properties of two types of the effective potentials. In 
Sec. Ill, the CMD method is briefly surveyed. We then newly introduce the EPAC method and present its numerical 
implementation. In Sec. IV, the results of the numerical tests for the accuracy of the EPAC method are shown in 
comparison with the CMD and the exact results. A possible improvement of the EPAC method is also discussed. The 
conclusions are given in Sec. V. 

Throughout this work, we consider a quantum particle of mass m in a one-dimensional potential V{q) at temperature 

T. 
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II. DEFINITIONS OF THE EFFECTIVE POTENTIALS 



A. Effective classical potential 

At first, we summarize Feynman's definition of the effective classical potential |f|. The quantum canonical 
partition function for a one-dimensional system at temperature T is expressed in terms of the imaginary time path 
integral 



Z a = 



30 j.q(ph)=q 

dq / Vq exp 

oo J q(Q)—q 



where (3 = 1/kgT and Sslq] is the Euclidean action functional 



~m<f + V(q) 



SE[q] = dr 
Jo 

The imaginary time average value of q{r) is defined as 

1 f (m 

qo = mJo dTq{T) - 



(i) 



(2) 



(3) 



This is the zero mode in the Fourier modes of q{r), q n = (l/(3h) J^ h d Te t ( 2 ^ n /i 3h )' r g(r), and is referred to as the path 
centroid. Inserting 1 = dq c 8(qo — q c ) to the integral of Eq. JQ), we obtain 



Z 8 = 



f q(f3h)=q 

dq c I dq Vq 5(q Q - q c ) e~ Se/H 

-OO J 9(0)— q 



dq c pp(qc) 



Here the centroid density 



dq / Vq S(q - q c ) e~ Se/ \ 

— 00 J q(0)—q 



and the effective classical potential 



log 



2n[3h 2 



P%{lc) 



(4) 



(5) 



(6) 



have been introduced as the functions of position centroid variable q c . This type of effective potential, Eq. 10, contains 
the effects of quantum fluctuations of q n ^o modes and is also called effective centroid potential 28] , constraint effective 
potential 22], or Wilsonian effective potential 
the PIMC or the PIMD techniques |2ll EE E3 



37 , |38| • This effective potential can be numerically evaluated by using 
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B. Standard effective potential 



In this subsection the effective action formalism is briefly reviewed and the standard effective potential is then 
intro duced mean s of the derivative expansion of the effective action. All the contents shown below can be seen in 
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Consider an imaginary time quantum theory in the presence of an external source J(t). The quantum canonical 
partition function of a one-dimensional quantum system with J(t) is expressed as 



2/3 [J] = 



dq 



q(0h)=q 



T>q exp 



'-oo J q(o)=q 
while the generating functional is defined as 



1 



~ h SE[q] + h ' d T J { T ) ( l{ T ) 



W P [J\ = %]ogZp[J\. 



(7) 



(8) 



The functional derivatives of W|g[J] with respect to J(t) produce the connected Green functions in the presence of 
J(r) 



(9) 
(10) 



where T represents the time-ordered product. Then the imaginary time connected Green function G™ nn (r — s) in 
thermal equilibrium is obtained when the external source J vanishes 



S 2 Wp[J] 



Gp(r-a)-(4{r))p(q{a))p, 



,7=0 



where Gp(r — s) is the imaginary time Green function (the Matsubara Green function) 

G (r-s) = (Tq(T)q(s))f3 = 9(T-s) (q(r)q{a))p + 6(s - r) (q(s)q{r))p. 



(11) 



(12) 



Note that the expectation value (q(r))p is independent of the imaginary time in thermal equilibrium. The effective 
action T l p[Q] is defined by the Legendre transform of Wjg[J] 

r 0h 



This satisfies the equation 



Tp[Q] = -W [J] 
8Tp[Q] 



cItJ(t)Q(t) 



(13) 



(14) 



SQ(t) 

which becomes the quantum mechanical Euler-Lagrange equation in the J — > limit, i.e., the principle of the least 
action. The functional derivative of Eq. I|14f> with respect to J(s) leads to 



I 



<* PW P [J] PT P [Q] = 

5J{s)6J(u) SQ{u)6Q(t) 1 



When the Fourier mode expansions with the Matsubara frequencies ui n — 2irn/ ' (ih are defined, we get 

1 °° 

( r -- s ) = jfi }^ e g p ( w «)' 



n— — oo 
oo 



SQ(t)SQ(s) [3h n £^ x " 5QSQ 
Then a general relation is obtained from Eq. (|15fl in the J — > limit 



(15) 

(16) 
(17) 

(18) 
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Thus the imaginary time (two-point) Green function is obtained with the knowledge of the effective action T^Q]. 
In a similar way, the imaginary time n-point Green function can be obtained with the knowledge of nth functional 
derivative of the effective action. Therefore, the effective action formalism reviewed here provides a powerful scheme 
to calculate the correlation functions. 

In general, the effective action rg[Q] is nonlocal in the imaginary time r, so that it is formally possible to expand 
r / g[Q] in a series of the terms involving imaginary time derivatives of q(r) (the derivative expansion) 

r(3K r i 

Vp(Q) + -Z P {Q)Q 2 + 0(d 4 ) 



rpn. 

I>[Q] = / dr 
Jo 



2 



(19) 



where Vp{Q) is called the standard effective potential, the leading order of the derivative expansion of I^gfQ]. 

A significant feature of Va(Q) is its convexity, which can be easily shown below. When we set Q{t) — > Q: constant, 
the effective action is written as I^fQ] = phVp(Q). In a similar way, when we set J(t) — > J: constant, the generating 
functional Wjg[J] can be written in terms of the density w(J): Wp[J] — j3Hw(J). Then Eqs. Ijl4|l and (|15|l become 

dV ^ = J, (20) 



dQ 

1 d 2 w p {J) d 2 V p {Q) 



(3% dJ 2 dQ 2 



= 1. (21) 



Considering the inequality 



one obtains 



— 0j2 — - ^ { (go - {qo)p) )p > U) (22) 



Therefore, the standard effective potential Vp{Q) is always convex, while the effective classical potential Vp(q c ) is not 
necessarily convex. 

C. Relationship between Vp(q c ) and Vp(Q) 

In this subsection we mention the relationship between the effective classical potential Vp{q c ) and the standard 
effective potential Vp(Q) EE1 EE Q ■ Replacing 5(q — q c ) in Eq. J5J) by its integral representation 

o pc+ico 

S(qo~q c ) = i- / dJe^*-««), (24) 



2m 



one obtains 



lOirRtfl r r c+io ° r°° ri( l3h ^ = i 



m 27ri 



c— -ioo J — og J q(0)— q 



m 27TZ 



(25) 



c— zoo 



When the low temperature limit (3 — ► oo is taken, the integral J d J can be evaluated by the saddle-point method 

e -PV g [q c ) = q e ic = -££ 

= Ce"^ w , (26) 
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where C is a constant. Therefore these two effective potentials are equal 

lim V p c (q c ) = lim V (q c ), (27) 

except for an additive constant term. Namely, both the effective potentials coincide with each other at the zero 
temperature. This relation also ensures the convexity of the effective classical potential VS (q c ) in the low temperature 
limit. 



III. REAL TIME CORRELATION FUNCTIONS 

Now we proceed to the calculations of the real time quantum correlation functions starting from the effective 
potentials. First we briefly review the CMD method [3] which is an approximation using the effective classical 
potential Vp(q c ). Then we newly introduce our EPAC method, which is a novel approximation using the standard 
effective potential Va(<5)- In this section, we concentrate on the two-point position correlation function 

Cp{t) - {q(t)Q(0))p, ( 28 ) 

to clarify the essence of the approximations. 



A. Centroid molecular dynamics method 

In the CMD, a real time classical equation of motion for the position centroid variable q c (t) on the effective 
classical potential V^(q c ) is introduced 



;q c (t) = FS(q c ) = - 



dVp C (qc 
dq c 



(29) 



The centroid force FZ(q c ) is evaluated by the Ehrenfest relation for q c Q 



13 

F${q c 



1 



dq / Vq S(q - q c 

Pp J-oo Jq(Q)=q 



(3h J Q dq(s) 



■ S E /h 



_ I l_ f h dV(q) Y 
" \0hj o dq{s)] 



(30) 



where 



denotes the quantum mechanical average with a constraint qo — q c . Using the centroid trajectory q c (t) 



generated from Eq. I|29|) , one can construct the centroid correlation function in a classical fashion 

I 



C7|(t) = (q c (t)q c (0)) Pfi =2-1 d( Lc 3»(t)?c(0) Ppilc) 



1 



OO pOC 



oo <J — oo 



'P 

dq c dp ( 
2irh 



q c (t)q c (0) e -/3[P 2 c/2m+^(«e)]. 



(31) 



Here we introduced the momentum centroid variable p c (t). Cao and Voth proposed the CMD approximation |28(, in 
which an approximate relation holds between two correlation functions 



C£(t) * C° AN (t), 



(32) 



where CS AN (t) is the canonical correlation function appearing in the linear response theory j^] 



c ? (t) = p J dxm-mx)q(o))p- 



(33) 
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In fact, the canonical correlation function is related to the real time quantum correlation function in the Fourier space 

Cp(u) = E(w) C$ an (oj), (34) 

where we introduced a function 

E(u) = ^(coth^ + l). (35) 

Therefore, the CMD method enables us to evaluate quantum dynamics from the correlation function Cg(i) calculated 
in a classical manner. The relation Eq. (|32|l is exact for a harmonic oscillator. For general potentials, Cp is exactly 
equal to C^ AN in the classical limit or for linear operators at t = 0, as shown in Appendix A. 

For general correlation functions of nonlinear operators, e.g. (q n (t)q n (Q)) p : the relation Eq. ilMl'l) is not valid even 
at t = 0. Reichman et al. have investigated this problem to propose the use of the higher-order Kubo transforms [sflj . 



B. Effective potential analytic continuation method 

In this subsection we newly present our quantum dynamics method. We begin with the review of the standard 
procedure of the analytic continuation from the imaginary time to the real time . Let us start from the Fourier 

coefficient of the Matsubara Green function [Eq. I|12l) ] 

,/jft 

Gp(u n ) = / dre^Gpir). (36) 
Jo 

Next the real time quantities, the retarded and advanced Green functions, are introduced 

cf(t) = e{t) < mum (37) 

Gj{t) =6{-t) ( [g(0),g(i)] h, (38) 

where [ , ] denotes the commutator. The Fourier coefficients of these Green functions are Gp(w) = dt e lult G^{t) 
and Gp(w) = f^dt e luJt Gp(t), respectively. One also introduces a couple of correlation functions, C^(t) = 
C p (t) = (q(t)q(0)) fJ , C<(t) = {q(0)q(t)) fj , and their Fourier coefficients, C>(lu) = J^dt e luJt C>(t), Cf{uo) = 
dt e ZUJt Cp(t). When one defines the spectral function 

/oc 
dte iut {[q(tU{0)])p (39) 
-OO 

= c>H-c<H (40) 

= (l-e-^)C>M, (41) 
then the quantities G^(u n ), G^ (u;), and G^(co) are expressed in terms of it 

/OO J / -t 
%r T——7 Pf>(<*')> ( 42 ) 

/OO J I -i 
.^^7TT e »^ < 43 > 

/OO J I -1 
IT — PpW), (44) 

where e is a positive infinitesimal. Therefore one can obtain the real time quantities G^{lo) and Gp(uj) from the 
imaginary time quantity Gp{u> n ) by means of the analytic continuation: 

Gf(u) = -i G^k)|^ , Gf(u) = -i G (u n )\ iu)n ^_ ie . (45) 
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Since the spectral function pp{to) is also expressed as 

and the real time quantum correlation function Cp(t) = (q(t)q(0))p is expressed as 



2tt 



2vr 



1 



1 



(46) 



(47) 



one can obtain the real time quantum dynamics with the knowledge of the imaginary time quantity via the analytic 
continuation procedure. 

Now we present the EPAC method, with which one can readily perform the analytic continuation procedure shown 
above, utilizing the effective action formalism [3,13 EH 13- As seen m Eq. 1|18|) . the imaginary time two-point 
connected Green function is expressed in terms of the effective action r / g[Q]. Now we employ the local potential 
approximation 



i-pn 

Tf}[Q] = / dr 
Jo 



1 



S2 



Vp(Q) + -mQ 



(48) 



where any derivative terms are dropped except for the kinetic term. From Eq. (|48|l . the second functional derivative 
of rg[Q] becomes the second derivative of the standard effective potential Vp(Q) 



SQ(t)5Q(s) 



= —m- 



dT 2 



dQ 2 



S(t-s), 



-(uj n ) = mu) n + 



(49) 
(50) 



SQSQ v ~" / '""~ n ' dQ 2 ' 

Note that d 2 Vp/dQ 2 must be evaluated at Q = Q m i n which gives the minimum of Vp(Q). This is because Eq. I|18|l 
holds in the J — » limit and the value of Q is fixed at Q = Q m i n in Eq. Q20|l . ft should be also noted that Q m in is 
independent of the imaginary time in thermal equilibrium. When we define the effective frequency 



CO/3 



1 d 2 Vp 



m dQ 2 



Q=Q mi „ 

from Eqs. Ijllfl . (|18|l . (|5U|I . and (|51|l . the Fourier transformed Matsubara Green function G^(u„) is written as 

h 



mu 2 + mui 



(51) 



(52) 







Here wc omitted the Fourier coefficient of the constant term (5^ in because it has nothing to do with the procedure 
described below. 

Following Eq. (|45ll . we obtain the real time quantities from the imaginary time quantity Gp(u n ) by the analytic 
continuation: 

h 



h 



l 



l 



u) — top + it u> + up + it 
1 1 



Imujp 

Then, using Eq. the spectral function has the form 

nft 



LJ — LUp — it LJ + LOp — it 



mujp 



S(uj — cop) — S(lu + LOp) 



(53) 
(54) 

(55) 
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Using Eq. Q47JI. we obtain an approximate real time quantum correlation function Cp C (t) with the constant term 
O 2 ■ 

^5 min 



°/3 



(*) 



( h 



V 2muj 







■ coth j cos ujpt — i 



2 

min - 



(56) 



We call this novel procedure the effective potential analytic continuation (EPAC) method, an approximation based on 
the effective action formalism and the derivative expansion. The EPAC correlation function Cp C (t) is similar to the 
exact quantum correlation function of a harmonic oscillator. However, note that the frequency u)p is the result coming 
from the curvature of the standard effective potential Vp(Q) at Q = Q m i n - It should be further noted that the EPAC 
differs from the analytically continued effective harmonic theory |28j which is based on the centroid variables. In this 
effective harmonic theory, the optimized frequency ui c {q c ) appearing in the time correlation functions comes from the 
curvature of the effective classical potential V£(q r ). Therefore, such effective frequency u> c (Qc) can be imaginary in 
some cases because Vp(q c ) is not always convex |41|. 



imaginary time 
path integral 



PIMC/PIMD 



effective 
classical 
potential 


Legendre 
*~ 


standard 
effective 
potential 

V (Q) 


transform 



classical 
eq. of motion 



CMD 

-tCAN / 



analytic 
continuation 



EPAC 



real time quantum 
correlation function 

c P (t) = (mm)? 



FIG. 1: Schematic diagram of the calculations of the real time quantum correlation function from the effective potentials: The 
Cao-Voth CMD method and the presently proposed EPAC method. 



To perform the EPAC procedure, one needs to calculate the standard effective potential Vp(Q). Once the effective 
classical potential VS(q c ) is obtained, it is straightforward to calculate the standard effective potential as shown below 
El- 

From Eqs. (@J and JSJ, we have 



2nf3h 2 



d?ce /3[^.-v#(g.)] j 



with a constant source J. Therefore, we obtain wp{ J) from V^(q c ), 



dq c e 



c. 



(57) 



(58) 
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where C is a constant. As in Eq. I|13(l . Vp(Q) is obtained through the Legendre transform 

V P {Q) = sup {JQ - wp(J) }. (59) 
j 

Here we find the maximum value of JQ — Wfj{J) with varying J for each given value of Q. This definition is equivalent 
to Eq. (|13fl though this is more suited for numerical evaluation. 

In Fig. ^ we show the flow chart of the calculations of the real time quantum correlation function by the CMD 
method and by the EPAC method. 

Finally we comment on the other contributed works of quantum dynamics methods based on the effective potentials. 
One of them is the direct use of the equation of motion Eq. I|14|) in the J — > limit 0[. Although this includes 
difficulty that Eq. (|14fl contains (infinitely) many time derivatives, more detailed analyses have been performed |43| . 
On the other hand, it has been known that, in the zero temperature limit, the CMD equation of motion [Eq. 
can be obtained from Eq. (|14J) in the J — > limit, employing the leading order derivative expansion and the analytic 
continuation. The meaning of the CMD equation in the zero temperature limit has been discussed in connection with 
the minimum energy wave packet [44I , while other related topics have also been discussed in Ref. 

IV. NUMERICAL TESTS 

In this section we test the accuracy of the EPAC method in a one-dimensional system by calculating the real time 
quantum correlation function C^ c (t) [Eq. I)5(i|l] . We compare C^ c (t) with the exact real time quantum correlation 
function 

c &) = ^T l J2 e - 0Ene ' i{Em - En)t/h \( m ^ n )\ 2 - (so) 

^ n m 

For comparison, we also show the CMD results of the centroid correlation function Cp(t) [Eq. Q31JI]- In this case, 
C C (t) is compared with the exact canonical correlation function 



1 _ e -(3(E m -E n ) 



e 



< E ™- En W h \(m\q\n)\ 2 . (61) 



The exact results of Cp(t) and C^ AN (i) are obtained from the eigenstates and the eigenvalues by solving the 
Schrodinger equation numerically. Throughout the numerical evaluation, we employ natural units, m = h = ks = 1. 
We consider a quantum particle moving on a double well potential 

V{q) = -\q 2 + ^q\ (62) 

All the calculation have been performed at (3 = 1 and 10. This double well potential model has already been 
analyzed in detail by use of the CMD and the other related methods 0|. The parameters chosen here correspond 
to, for example, a proton potential where the potential height AV is about 0.5 kcal/mol and the length between the 
potential minima Aq is about 1 A; this is a relatively shallow double well potential model. The temperatures in the 
present study are T = 1.6AV ~ 400[K] for j3 = 1 and T = 0.16AV ~ 40 [K] for = 10. 

In the following subsections, first we show the CMD results with the exact canonical correlation function. Next the 
EPAC results are compared with the exact quantum correlation function. We also present the Fourier transformed 
correlation functions to indicate the difference between the CMD method and the EPAC method clearly. 

A. Centroid molecular dynamics method 

For the CMD calculation, at first we need to evaluate the centroid force F^(q c ) [Eq. (J30J] before solving Eq. i(2*9"|) 
numerically. For this purpose, we employed the PIMD technique. The basis of the PIMD technique is the discretized 
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FIG. 2: /3-dependence of the effective classical potential Vp(q c ). We set V^(0) = 0. 
representation of the quantum canonical partition function [Eq. l(T|l] 



(64) 



where P is the Trotter number and 3>p(q) is the potential of quasi-particles q = {q±, .., qp}, 



*p(q) = E 



v mP , „ 1 



with the periodic boundary condition gp+i = gi. By performing the constant-temperature molecular dynamics (MD) 
for the quasi-particles, static properties corresponding to the quantum canonical partition function Zp can be obtained 
as the averages over such MD-based configurations [2j . In this framework, the centroid force 





Qj - 1c) 




p 


dV(q 3 )~ 
1 dqj 


e -/3*i»(q) 


I---IUUd qj 6 




Qj - 9c) 


e -/3* P (q) 



^ C (9c) = - p hm ^ — ^ _ D X 1 ~ ^ , (65) 



was computed at each fixed centroid position q c by use of the normal mode PIMD algorithm |28|, |47j including the 
Nose- Hoover chain (NHC) thermostats 0], which ensure the generation of the static canonical distribution. We used 
10 7 configurations to evaluate F^(q c ) for each of 51 grid points of q c , equally spaced between [—2.5,2.5]. Then we 
fitted the centroid force F^(q c ) to the 25th polynomial function. After this fitting, the effective classical potential 
Vp(q c ) was obtained from the integration of FS(q c ). Figure El shows the classical (bare) double well potential V{q) 
and the effective classical potential Vp(q c ) at (3 = 0.1, 1, 10 and 100. As shown in Sec. II [Eqs. lO and VS{q c ) 
is not convex at higher temperature, though it becomes convex as the temperature lowers. 

The second step is to solve the CMD equation of motion [Eq. (|29|) ] numerically. The centroid trajectories needed to 
compute the centroid correlation function C%(t) were calculated by means of the dual sampling procedure j^EUIS- 
In this procedure, at first the sampling of the initial position and momentum centroids was performed using the CMD 
with the centroid-attached NHC thermostat. And then the centroid dynamics Eq. I|29(l was computed from the initial 
distribution in a microcanonical manner to obtain the centroid trajectories. The centroid correlation function was 
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FIG. 3: The exact canonical correlation function C^ AN (t) and the CMD correlation function Cp{t): (a) for /3 = 1 and (b) for 
P = 10. 



evaluated from such centroid trajectories spanning 10 6 steps. This is the correct evaluation of Cp(t) following its 
definition Eq. (f?TT|) . 

Figure shows the exact canonical correlation function C^ AN (i) and the centroid correlation functions C|(i) at 
two temperatures, (3 = 1 and 10. In both cases, we can see that the CMD approximation well reproduces the exact 
short time behavior, as was already found in Refs. and El This is because the CMD approximation Eq. Ij^2(l is 
exact at t = (see Appendix A). However, at t > the relation Eq. (|32jl does not exactly hold, so that the long time 
CMD behavior eventually deviates from the exact one, even though the centroid force F^{q c ) is calculated correctly. 
For (3 = 10, the coincidence of C^(t) with the exact C9 AN (i) persists up to t ~ 4, whereas it breaks down at t ~ 2.5 
for (3 = 1. Namely, the CMD result for /3 = 10 is better than for (3 = 1. This is because, for such shallow double well 
potential as we used here, the shape of V^(q c ) changes into quasi-harmonic one as the temperature lowers (see Fig. 
0); this results in a weaker ensemble dephasing effect |30j,|32|. However, this property does not hold for deeper double 
well potentials whose effective potentials have strongly anharmonic shape |43, |50j . 



B. Effective potential analytic continuation method 



Following the procedure presented in Sec. Ill B, we have calculated the standard effective potential Vp(Q) from 
the effective classical potential Vg (q c ) which was obtained from the PIMD calculations in the last subsection. Figure 
0] shows Vp(Q) at (3 = 0.1, 1, 10 and 100. All of the curves are convex, as analytically shown in Eq. (j23(l . Especially 
for (3 = 100, Vp(Q) has almost same potential shape as Vp(q c ) in Fig. as expected from the relation Eq. 127(1 . 
As for the high temperature limit, V^(q c ) and V^(Q) have quite different shape. For (3 — > 0, the effective classical 
potential Vp(q c ) is always equal to the classical (bare) potential, because there is no quantum correction to V^(q c ) 
in the high temperature limit. On the other hand, the standard effective potential Vp{Q) generally changes from the 
classical one even in the high temperature limit, because Vp(Q) receives the corrections from the quantum fluctuation 
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FIG. 4: /3-dependence of the standard effective potential Vp(Q). We set 1^(0) = 0. 

and the thermal fluctuation. In fact, for the double well potential considered here, as suggested in Fig. QJ Vp(Q) in 
the [3 — > limit becomes a singular potential which has an infinite effective frequency ujp — > oo. As an exception, it is 
well known that for the harmonic oscillator, Vp(q c ) and Vp(Q) have the same potential shape at all the temperatures, 
because in this case the quantum and/or thermal fluctuation contributes only to the corrections to the constant terms 
in Vp(q c ) and Vp(Q). It should also be noted that Q m in which gives the minimum of Vp(Q) is always zero for the 
classical potential Eq. (ffi2*)l considered here, because the potential Eq. (ffi2*)l is Z 2 (q <-> —q) symmetric and Vp(Q) 
always has a convex (d 2 Vp/dQ 2 > 0) shape |51| . 

In Fig. we show the real part of the exact quantum correlation function Cp(t) and the real part of the calculated 
EPAC correlation function Cp C (t) at (3 = 1 and 10. At the higher temperature [Fig. Eta)], the EPAC method 
reproduces well the exact correlation function at short time less than t ~ 2; such a good agreement in the short time 
behavior is similar to the results of the CMD method. However, in contrast to the CMD, the approximate function 
C^ c (t) deviates from the exact quasi-periodic behavior at longer time because of its completely periodic oscillation. 
This is caused by the fact that Cp G {t) has no dephasing because it consists of only one oscillation mode with the 
frequency up [see Eq. (|56|) ]. On the other hand, for j3 = 10, the initial value Cg C (0) slightly deviates from the exact 
one, while Cp C (t) reproduces the global oscillating behavior of Cp(t) better than the CMD correlation functions. 
This initial value deviation is because the EPAC approximation [Eq. I|56|l] is not necessarily exact at t = 0. The 
initial value C^ c (0) can be improved by employing the higher order derivative expansion (see the next subsection); 
it should become exact in an infinite order expansion. The EPAC is therefore a method to capture the oscillating 
behavior of quantum systems; it should be very effective for systems in which quantum coherence is significant. 

The essential features of the correlation functions can be seen more clearly in the Fourier space. The Fourier 
component of the exact real time quantum correlation function is 

/OO 1 
e^Cp(t) = — V V e-^-(27r) *(w - w m , n ) \(m\q\n}\ 2 , (66) 

where uj n = E n /h and uj m ^ n — (E rn — E n )/h. The correlation function Cp{u>) has a discrete spectrum and each oj m ^ n 
mode contributes to Cp(u>) with the weight e~^ hu)n |(m|g|n)| 2 . Therefore, for finite /3, Cp(t) exhibits quasi-periodic 
oscillation as a result of the interference between many u> m ,n modes. As for the canonical correlation function, the 
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FIG. 5: The real part of the exact quantum correlation function Cp(t) and the real part of the EPAC correlation function 
C0 C (t): (a) for = 1 and (b) for = 10. 

Fourier component is Cg AN (w) = Cpluj) / 'E(u>), as expressed in Eq. I|34l) . Since the function E(u>) is unity at uj = 
and is a monotonously increasing function of ui, the higher frequency modes in Cp AN (ui) are less intensive than those 
in Cp(u). However, it is noted that the difference between Cp(u) and (7g AN (u;) lies just in intensity; the discrete 
spectrum coming from the delta-function in Cp(u>) [Eq. (|66|l ] is also the case of Cp AN (iv). 

On the other hand, for the CMD correlation function Cp(t), its Fourier spectrum is, in general, continuous because 
there is no origin of a discrete spectrum in the classical definition of Cp(t) in Eq. I|31fl . Therefore, the thermal 
fluctuation causes the destructive interference between an infinite number of oscillating modes; this is the ensemble 
dephasing. This effect leads to the result that the CMD method underestimates the exact quasi-periodic oscillation. 
Only the harmonic oscillator is free from this dephasing due to the interference. As for the EPAC correlation function, 
however, its Fourier component has the form 

/oo 
e lu)t C^ c {t) 
-OO 

= E{u) -V- [ S(cj + cjp) + S(oj - up) } + (2tt) S(u)Q 2 min , (67) 

which consists of a single oscillation mode with the frequency ujp. Evidently, Eq. (|67J) gives a discrete spectrum as 
well as Cp{uj). 

Figure |B| shows the Fourier transformed correlation functions, (7g AN (a;) and C|(o;). The spectra are plotted after 
reducing artificial Fourier ripple due to the finite time cutoff. The exact poles at u> = uj m ^ n and the EPAC pole at 
u; — cup are plotted in the same figure. As for the Fourier transformed canonical correlation function C9 (u>), at 
both temperatures, the peaks exist at the locations of a finite number of poles u m ,n- Hence C^ AN (t) is dominated 
by a limited number of oscillation modes, leading to the quasi-periodic oscillatory behavior. On the other hand, the 
CMD correlation function Cp{u>) has continuous spectrum for each temperature, while the EPAC correlation function 
has a discrete spectrum with a peak at oj — up. The CMD and the EPAC methods thus make a remarkable contrast 
with each other. For lower temperature (3 = 10, in Fig. HJb), fewer oscillation modes are dominant in C9 AN (o;) and 
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the EPAC pole ojp exists near the exact first pole 0*1,0, while the CMD correlation function C%(u>) still has a broad 
spectrum. 

As mentioned above, the EPAC is an approximation to extract the effective periodic oscillation from the exact 
correlation function. Therefore, it can be useful for the investigation of the systems where quantum coherence is 
significant. On the other hand, the CMD is a method to approximate the exact spectrum C^ AN (w) by a continuous 
spectrum, so that it is not proper to analyze the long time behavior of such systems. However, it has been found by 
Krilov and Berne, that the CMD gives much better results in a system under isobaric conditions, where the exact 
correlation function has a continuous spectrum [3^]. They also suggested that the CMD might give fairly accurate 
results in a dissipative system, because in such a system the quantum coherence is dcphased through the interaction 
with the dissipative environment not to play an important role at longer time. 



C. Possible improvement 

The EPAC method can be improved by means of the higher order derivative expansion of the effective action 
r^[Q], For example, the second order derivative expansion in Eq. (flTTf) is 



T fi [Q] = / dr 
Jo 



V P (Q) + -Z fj (Q)Q 2 



(68) 



Note the difference from Eq. (|48|) . Then we obtain the correlation function as a result of the second-order EPAC 
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FIG. 7: The real parts of the quantum correlation functions at zero temperature: the exact quantum correlation function C(t), 
the leading order EPAC correlation function C AC (t), and the second-order EPAC correlation function C AC ^(t). 



where Zp = Zp(Q min ) and 



1 d 2 Vp 



Zp dQ 2 



(70) 



The second-order EPAC correlation function again consists of a single oscillation mode, because the second-order 
derivative expansion only introduces the quantum/thermal correction to the particle mass, m — > Zp, and never 
changes the number of the poles. In fact, it is expected that the n(> 2)th order derivative expansion improves the 
qualitative behavior of the EPAC correlation functions by introducing additional poles. In such cases, the nth order 
EPAC correlation function Cp (t) would exhibit quasi-periodic oscillation more similar to the exact behavior. 

To see the effectiveness of this improvement, we now present a zero temperature example. At zero temperature, 
the exact real time quantum correlation function is 



C(t) = <0|g(t)g(0)|0> 



< E *»- E °W h \(m\q\0)f 



(71) 



On the other hand, the leading order EPAC correlation function is 

h 



Q 



2 

min J 



where uj e g 



. 1 - 

wp\p^oa an d the second order EPAC correlation function obtained from Eq. I|69|l is 

n 



C AC(2) W 



2Z e gLO eS 



^ + Q 



2 

min 7 



(72) 



(73) 



and wffr 



uiplp^oc- To evaluate these zero temperature quantities, ui c s, Z c g, and u>^ s , we 



where Z eS = Zp\p_ >QO lu cS 
employed the renormalization group technique |37j for convenience; the effective frequency ta e g was calculated by 
solving the local potential approximated Wegner- Houghton equation [ssi 1^1 . while for Z c g and uj^ s , we used the 
values obtained by solving the proper time renormalization group equation |50J . Figure [7| shows the real part of the 
exact quantum correlation function C(t) and the real parts of the EPAC correlation functions: the leading order 
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C AC (t) and the second order C AC{2 \t). It is seen that the deviation of the leading order EPAC correlation function 
C AC (t) from the exact quantum correlation function is improved in the second-order EPAC (C AC ^ 2 \t)), both in its 
initial (t — 0) value and in its long time behavior. 

The higher order derivative expansion is required especially for the analyses of many-body systems at higher 
temperature, because the multipole contribution to the correlation functions is significant for describing the decay 
of them. However, it is not easy to compute the higher derivative terms in Tp, especially at finite temperature. 
Therefore, it is a future task to establish an algorithm for computing them. 



V. CONCLUDING REMARKS 

As a novel approximation method to evaluate the real time quantum correlation functions at finite temperature, 
we have newly proposed the EPAC method. The EPAC method has been tested in the one-dimensional symmetric 
double-well system, in comparison with another approximation scheme, the CMD method. The EPAC method and 
the CMD method are based on each different type of effective potential, the standard effective potential Vp(Q) and 
the effective classical potential VS(q c ), respectively. At first, we have evaluated Vp{q c ) by means of the normal mode 
PIMD calculation, and then Vp(Q) have been obtained from it. It has been found that these effective potentials have 
the predicted formal properties: The equivalence of them in the zero temperature limit and the convexity of Vp{Q). 
Then the real time two-point position correlation functions have been calculated by means of each approximation 
method. The CMD approximation is found to be good at short time range owing to its exactness at t = 0, while at 
longer time it largely deviates from the exact correlation function because of the ensemble dephasing. On the other 
hand, our EPAC approximation can reproduce the long time oscillating behavior which originates from the quantum 
coherence of the system. Therefore, the EPAC should be very effective for the system in which quantum coherence is 
significant, such as the quantum double well system considered here. We have also suggested that the EPAC method 
can be improved by the higher order derivative expansion, and have shown, as an example, the result of the second 
order improvement at zero temperature. It has been seen that this improvement procedure works very well for this 
example. 

In this paper we have restricted the arguments to the evaluation of the real time two-point position correlation 
function Cpit) = (q(t)q(Q))p. Also for general correlation functions of nonlinear operators, e.g., (q n (t)q n (0)}p, the 
EPAC can be applied to evaluate them, because the re-point correlation function is obtained with the knowledge of 
the nth functional derivative of the effective action [s^l li^l Isfil Isfif . This will be an interesting subject to be examined 
in near future. 

As matters of theoretical interest, there are a number of subjects to be investigated. The role of the standard 
effective potential in the quantum transition-state theory should be clarified. Furthermore, it is interesting to 

find possible relationship between the effective potential-based methods of the present type and the other approaches 
to quantum dynamical correlations |53|. It is worthwhile to pay attention to the fact that, for example, the effective 
potential-based quantum dynamics has been discussed in the context of the projection operator approach \5^ . 

Finally we mention the direction of the applications of the EPAC to real molecular systems of chemical interest. 
From the results given in Sec. IV, it is suggested that the EPAC should be very effective for the systems where 
quantum coherence plays an important role. For example, proton transfer reactions [fifij are typical phenomena of 
such category, because the small mass of proton makes quantum coherence significant even at room temperature. The 
EPAC will properly represent the quantum oscillating behavior accompanied with proton transfers. For example, for 
a molecular reaction system where the intrinsic reaction coordinate (IRC) [3] is well-defined, once the IRC and the 
potential energy surface along it are provided, it is straightforward to calculate the approximate real time quantum 
correlation function by use of the EPAC method. However, for reactions occurring in solvents, not only the IRC 
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but the full quantum calculations treating many degrees of freedom should be implemented. For the application of 
the EPAC to such many-body systems including quantum coherence, more efficient sampling algorithms or novel 
approximation schemes must be developed to calculate the effective classical potentials needed in the EPAC. 
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APPENDIX A: ON THE RELATION BETWEEN CJ(0) AND C^ AN (0) 

It has been known that the CMD approximation [Eq. is exact at t = 0, i.e., C|(0) = C^ AN (0) 2 

Note that this relation holds only for the correlation function of linear operator, as seen in the following simple 
derivation: 
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and the definition of path centroid go = q(^>o)/ /3fo- 
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